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Abstract: We consider the resummation of soft gluon emission for squark and gluino 
hadroproduction at next-to-leading-logarithmic (NLL) accuracy in the framework of the 
minimal supersymmetric standard model. We present analytical results for squark-squark 
and squark-gluino production and provide numerical predictions for all squark and gluino 
pair-production processes at the Tevatron and at the LHC. The size of the soft-gluon 
corrections and the reduction in the scale uncertainty are most significant for processes 
involving gluino production. At the LHC, where the sensitivity to squark and gluino 
masses ranges up to 3 TeV, the corrections due to NLL resummation over and above the 
NLO predictions can be as high as 35% in the case of gluino-pair production, whereas 
at the Tevatron, the NLL corrections are close to 40% for squark-gluino final states with 
sparticle masses around 500 GeV. 
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1. Introduction 

The search for supersymmetry (SUSY) [1, 2] is among the most important tasks at cur- 
rent and future colliders. Squarks and gluinos, the coloured supersymmetric particles, are 
expected to be produced most copiously in hadronic collisions. Searches at the proton- 
antiproton collider Tevatron with a centre-of-mass energy of y/S = 1.96 TeV have placed 
lower limits on squark and gluino masses in the range of 300-400 GeV [3, 4]. The proton- 
proton collider LHC with = 14 TeV design energy will extend the range of sensitivity 
to squarks and gluinos with masses up to about 3 TeV [5, 6, 7]. 
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In the minimal supersymmetric extension of the Standard Model (MSSM) [8, 9] with 
R-parity conservation, squarks and gluinos are pair-produced in collisions of two hadrons 
hi and hi'- 

h\h 2 -> qq,qq,qg,£ig + X . (1.1) 

In Eq. (1.1) and throughout the rest of this paper we suppress the chiralities of the squarks 
q = (qL,qn) and do not explicitly state the charge-conjugated processes. We include 
squarks q of any flavour except for top squarks. The production of top squarks [10] has to 
be considered separately since the strong Yukawa coupling between top quarks, top squarks 
and Higgs fields gives rise to potentially large mixing effects and mass splitting [11]. 

Accurate theoretical predictions for inclusive cross sections are crucial to derive exclu- 
sion limits for squark and gluino masses [3, 4] and, in the case of discovery, can be used to 
determine sparticle masses [12] and properties [13]. The cross sections for the squark and 
gluino pair-production processes (1.1) are known at next-to-leading order (NLO) in SUSY- 
QCD [14, 15, 16]. Electroweak corrections to the 0(o%) tree-level production [17, 18, 19, 20] 
and the electroweak Born production channels of 0(aa s ) and 0(a 2 ) [21, 22] are signifi- 
cant for the pair production of SU(2)-doublet squarks q~i and at large invariant masses in 
general, but they are moderate for total cross sections summed over all squark species. 

The NLO SUSY-QCD corrections to squark and gluino hadroproduction reduce the 
renormalization- and factorization-scale dependence of the predictions. In general these 
corrections also significantly increase the cross section with respect to the Born predic- 
tions [23, 24, 25] if the renormalization and factorization scales are chosen close to the 
average mass of the pair-produced sparticles. A significant part of these large corrections 
can be attributed to the threshold region where the partonic centre-of-mass energy is close 
to the kinematic threshold for producing massive particles. In this region the NLO cor- 
rections are dominated by the contributions due to soft gluon emission off the coloured 
particles in the initial and final state and by the Coulomb corrections due to the exchange 
of gluons between the massive sparticles in the final state. The soft-gluon corrections can be 
taken into account to all orders in perturbation theory by means of threshold resummation. 

Previous work has addressed the soft-gluon resummation for squark-antisquark and 
gluino-gluino production at next-to-leading- logarithmic (NLL) accuracy [26, 27]. For the 
squark-antisquark production process the dominant contribution to the next-to-next-to- 
leading order (NNLO) correction coming from the resummed cross section at next-to- 
next-to-leading-logarithmic (NNLL) level has been studied in [28]. Moreover, a formalism 
allowing for the resummation of soft and Coulomb gluons in the production of coloured 
sparticles has been presented in [29, 30], and bound state effects have been studied for 
gluino-pair production in Ref. [31]. Additionally, threshold resummation for single colour- 
octet scalar production at the LHC has been investigated in [32] . 

In this work, we present the analytical components needed to perform NLL resumma- 
tion for squark-squark and squark-gluino pair-production. In addition, we provide numer- 
ical predictions for the entire set (1.1) of pair-production processes of coloured sparticles 
at the Tevatron and the LHC. 

The paper is structured as follows. In section 2 we review the formalism of soft-gluon 



- 2 - 



resummation. The calculation of the one-loop soft anomalous dimension matrices for the 
qq and qg production processes is discussed in section 3. We present numerical results 
for squark and gluino production at the Tevatron and the LHC in section 4 and conclude 
in section 5. A more detailed description of certain aspects of our calculation and some 
explicit formulae that enter the expressions for the resummed cross sections are collected 
in the appendices. 



2. Soft-gluon resummation 



In this section we review the formalism of threshold resummation for the production of 
a pair of coloured massive particles. Since the corresponding theoretical expressions have 
already been discussed in detail in Ref. [27], we shall be brief. 

The inclusive hadroproduction cross section a^^^ki for two massive SUSY particles 
k and I, where k, I can be a squark (q), antisquark (q) or gluino (g), can be written in terms 
of its partonic version (Tij->kl as 

<jh 1 h 2 ^ki{p,{ m2 }) = J dxidX2d ^{f~^^) 

x fi/hrixi, n 2 ) f j/h2 (x 2 , n 2 ) a i:j ^ k i(p,{m 2 }, fi 2 ) , (2.1) 



where {m 2 } denotes all masses entering the calculations, i,j are the initial parton flavours, 
fijh x and fj/h 2 the parton distribution functions, and p, is the common factorization and 
renormalization scale. The hadronic threshold for inclusive production of two final-state 
particles with masses 777,3 and m i corresponds to a hadronic center-of-mass energy squared 
that is equal to S = (7773 + 7774 ) 2 . Thus we define the threshold variable p, measuring the 
distance from threshold in terms of energy fraction, as 

(7773 + 7774) 2 



The partonic equivalent of this threshold variable is defined as p = p/(x\X2), where xi^ 
are the momentum fractions of the partons. This is a generalized version of the threshold 
variable used e.g. in Ref. [27]. It accounts for unequal masses of the pair-produced particles 
in the final state, making it applicable to the case of squark-gluino production. 

In the threshold region, the most dominant contributions to the higher-order QCD 
corrections due to soft gluon emission have the general form 1 

<log m /3 2 , 777 < 2n with j3 2 = 1 - p = 1 - ( m3 + m4 ) 2 ; (2.2) 

where s = X1X2S is the partonic center-of-mass energy squared. The resummation of the 
soft-gluon contributions is performed after taking a Mellin transform (indicated by a tilde) 



1 See section 3 for more discussion on the form of a threshold variable in the case of unequal masses. 
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of the cross section, 



Ch 1 h 2 ->kl 



(N, {m 2 }) = [ dp 



Jo 



P N 1 cr hlh . 2 ^ k i(p,{m 2 



}) 



(2.3) 



(N,{m 2 },p 2 ) . 



The logarithmically enhanced terms are then of the form a™ log m N, m < 2n, with the 
threshold limit (3 — > corresponding to N — > oo. The resummed cross section takes the 
schematic form [33, 34] 



in which all dependence on the large logarithm L = log N occurs in the exponent, and no 
term in the perturbative series P(a s ) grows with increasing N. Keeping only the g\ term 
constitutes the leading logarithmic (LL) approximation, including also the g 2 term is called 
the next-to-leading logarithmic (NLL) approximation, etc. Up to NLL accuracy it suffices 
to keep the lowest-order term in P. 

The all-order summation of such logarithmic terms depends on the near-threshold 
factorization of the cross sections into functions that each capture the effects of classes of 
radiation effects: hard, collinear (including soft-collinear), and wide-angle soft radiation 



where we have introduced the hard scale Q 2 = (7713 + 7714 ) 2 . Before we comment on 
each function separately, we recall that soft radiation is coherently sensitive to the colour 
structure of the hard process from which it is emitted [39, 35, 40, 36]. The various structures 
are labelled by the indices /, J in a way made more precise further below. 

The functions Aj and Aj sum the effects of the (soft-) collinear radiation from the 
incoming partons. They are process-independent and do not depend on the colour struc- 
tures. They contain the leading logarithmic dependence, as well as part of the subleading 
logarithmic behaviour, and are listed e.g. in Ref. [27]. 

The function Hij^kl,Jl incorporates only higher-order effects of hard, off-shell partons 
and therefore does not contain log N dependence. This hard function depends on the colour 
representations of the external particles in the partonic process. There are usually multiple 
tensors cj that can connect these colour representations, where / labels the possible tensors. 
For instance, in the case of squark-antisquark (with colour indices 03,04) production by 
the annihilation of light quarks (with colour indices 01,02) there are two colour tensors, 
which may be chosen as 



&hih 2 ^kl{N) = exp Lgi{a s L) + g 2 {a s L) + ... x P(a s ) , 



(2.4) 




l,Q 2 ,p 2 )A j (N + l,Q 2 ,p 2 ) 

(N, {m 2 },p 2 ) Sij^ k ijj(Q/ (Nfi),p 2 ) , (2.5) 



01(01,02503,04) — S aia2 S asa4 



(s- channel singlet), 



02(01,02503,04) = T£ ia2 TZ 



(s-channel octet) . 



(2.6) 
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The hard function jj is a matrix in this colour-tensor space, with the indices JI 

indicating the colour structure. Note that we paired the indices in example (2.6) according 
to the s-channel. Other choices are possible as well [39, 40, 37], but choosing an s-channel 
basis will be convenient at threshold. 

The soft function Sij^kijj in Eq. (2.5) is also a matrix in colour-tensor space, since soft 
emissions mix the connecting colour tensors. This soft function is constructed [40, 37] from 
an eikonal cross section, which in turn is defined in terms of the square of expectation values 
of products of Wilson-line operators belonging to the external particles in the process. 
These Wilson lines generate to all orders the soft-gluon radiation in the process and depend 
on the direction and colour representation of the corresponding external particle. To avoid 
double counting with the Aj and Aj factors in Eq. (2.5), the expectation values are divided 
by the square of expectation values of the Wilson lines themselves. In this way, collinear- 
soft radiation already included in the Aj and Aj factors is removed. What remains is a 
soft function whose perturbation series takes the form a™ log m iV, m < n, and therefore 
contributes only at NLL accuracy. 

Although the combination of the soft and collinear functions in the cross section is 
gauge invariant, the functions themselves are not automatically separately gauge invariant. 
The collinear functions only depend on the colour representations of the incoming partons. 
Therefore the gauge dependence of the soft function cannot depend on the colour structure 
of the process either. This implies that we can make the soft and collinear functions 
separately gauge invariant by rescaling them with a scalar in colour-tensor space. This 
rescaling has implicitly been performed in Eq. (2.5), where the soft function has been 

divided by yj S^ ng yj S^ lg as indicated by the bar on Sij^kijj- The factor S*" lg is the 
soft function for two incoming Wilson lines of flavour i and i annihilating into a colour- 
singlet 2 . By taking the square root of such a soft function, we effectively isolate the gauge 
dependence of a single line. Therefore this procedure works not only for qq or gg initial 
states but also for initial states that cannot annihilate into a colour-singlet, such as qg and 
qq. To compensate for the division factor in the soft function, the collinear functions Aj 



and Aj have been multiplied by the factors yj S^ n& and y S^ s respectively. Analytical 
expressions for these functions given in the literature (see e.g. Ref. [27]) explicitly include 
this multiplicative factor. 

Near threshold the soft function reduces considerably. For the inclusive cross section 
and our choice of colour basis, the matrix Sij^kijj becomes diagonal in colour-tensor space 
in the threshold limit (3 — > [27]. In this limit we have (suppressing particle flavour labels) 

Km Su(Q/(Nfi),fi 2 ) = SuSfjA^iQ/iNv),^) (2.7) 

with 

A<*>(C3/(iW) =ex P [/ 0/ "**M D ,l, (2 . 8) 



2 Note that if the colour representations are 3 and 3 this corresponds to the Drell-Yan process. For 
octets, it corresponds to Higgs production by gluon fusion. 
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where Sj°J is the lowest-order expression for the soft function, given by 



Sf] = tr(c\cj) . (2.9) 



The one-loop coefficients Dj are defined by 

Di = lim — 2Re(f//). (2.10) 

The values of the Dj coefficients for qq and qg production are calculated in section 3.4.3. 
The form of Eq. (2.8) follows from a renormalization-group equation for Sjj(Q/(Nfi)) 
[35, 37], with one-loop anomalous dimensions Tij^kljJ, often referred to as the "soft" 
anomalous-dimension matrix. If the calculations are performed in the axial-gauge with 
gauge vector n M , the one-loop anomalous dimensions are given by 

r« = r„ - J £ C 2 , p (l - log (2 { ^-f- ) - zvr) S U , (2.11) 

P={i,j} 

where the sum is over the two incoming particles, and \n\ 2 = —n 2 — ie, see Ref. [39]. The 
dimensionless vector v p is given by the momentum of the incoming massless particle p 
multiplied by y/2/s. The factors C^p are either Cf or Ca, depending on whether p is 
a quark or gluon, respectively. The subtraction exhibited in Eq. (2.11) results from the 

division by the factor described before. The matrix Tjj is the anomalous 

dimension matrix of the products of Wilson- line operators connected by the various possible 
colour tensors mentioned earlier. More details on its calculation are given in section 3.4. 
In the threshold limit the resummed partonic cross section becomes 

4%(N,{m 2 },^) = £ ^l klJ (N,{m 2 },^) C^ANAm 2 },^) (2.12) 



A t (N + 1, Q 2 , /x 2 ) Aj{N + 1, Q 2 , /i 2 ) A^ Hi/ (Q/(iv>), /i 2 ) , 



where cr^_^ kl j are the leading-order (LO) cross sections in Mellin-moment space. For the 
case of qq and qg production we present them in appendix A. The functions Cij^kij are 
of perturbative nature and contain information about hard contributions beyond leading 
order. This information is only relevant beyond NLL accuracy and therefore we keep 
Cij->kl,l = 1 in our calculations. 

Having constructed the NLL cross-section in the Mellin-moment space, the inverse 
Mellin transform has to be performed in order to recover the hadronic cross section 
u h\hi-^kl- in order to retain the information contained in the NLO cross sections [14, 15, 16], 
the NLO and NLL results are combined through a matching procedure that avoids double 
counting of the logarithmic terms in the following way: 

(NLL+NLO matched) / r 2i 2\ (NLO) / r 2i 2\ 

>>W },» ) = v y hlh2 S k i{p,{m },n ) (2.13) 



+ E / p~ n h /hl (N + i^ 2 )f 3/h2 (N + i^ 2 ) 

• • — J CT 



i,j=q,q,g 



~(res) /-.j r 2i 2\ ~ (res) / , r r 2t 2\ i 

(NLO) 



- 6 - 



We adopt the "minimal prescription" of Ref. [41] for the contour CT of the inverse Mellin 
transform in Eq. (2.13). In order to use standard parametrizations of parton distribution 
functions in x-space we employ the method introduced in Ref. [42]. 

3. Soft anomalous dimensions and Born cross sections for qq and qg pro- 
duction 

3.1 Kinematics 

To set the stage for the discussion of the soft anomalous dimensions we first introduce the 
relevant kinematical definitions that are used in the calculation. We consider the following 
generic process 

i(ai,Pi)j(a2,P2) -> k(a 3 ,p 3 )l(a4,p4) , (3.1) 

where the colour indices a« and the momenta of the particles pi are given in parentheses. In 
those cases where a final-state squark features in the process, summation over both squark 
chiralities (q L and q R ) and all possible squark flavours is implied, the latter being restricted 
by the choice of initial-state quark flavours. For the processes investigated here, i.e. squark- 
squark (kl = qq) and squark-gluino (kl = qg) production, top-squark final states are not 
possible since top quarks are excluded as initial-state partons. In view of the absence of 
top-squark final states, all squark-flavour and chirality states are considered to be mass 
degenerate with mass rriq. The gluino mass is denoted by trig. 

All analytical results presented in section 3 are derived for a general SU(iV c )-theory, 
with N c the number of colours. This means that the colour indices a« for gluons and gluinos 
can take N 2 — 1 different values, since these particles are in the adjoint representation. For 
(s)quarks, which are in the fundamental representation, the colour indices are -/V c -valued. 

The particle momenta featuring in the generic process (3.1) obey the on-shell conditions 
p\ = p\ = , p\ = m\ and p\ = m\. For the kinematical description of the reactions the 
standard Mandelstam invariants 

s = ( Pl + P2 ) 2 , t = (P1-P3) 2 and u = (pi - p 4 ) 2 (3.2) 

are used. In the centre-of-mass frame of the final-state particles the absolute value of the 
final-state momenta can then be written as 

\p3\cm = Ipilcm = ^k(3t/s, (3.3) 
with (3 = \J\ — (7713 + mi) 2 / s defined in Eq. (2.2) and 

K EE yi-^-^ 2 

The presence of the factor k is special to the case of unequal masses. As Eq. (3.3) 
shows, it occurs quite naturally in matrix-element expressions for the processes we consider 
in this study We could have defined the variable (3 1 = k(3 and taken moments with respect 
to this variable. Instead we have opted to use the variable (3 in our calculations in order to 
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(3.4) 



facilitate convolutions underlying the resummation. Because log/3' = log /3+log k, choosing 
f3' would have resulted in different subleading logarithmic terms. To NLL accuracy these 
differences in the expressions for the resummed partonic cross sections are cancelled by 
different terms arising from the convolutions. 

In order to present the results for the leading-order partonic cross sections it is helpful 
to introduce two more shorthand notations: 

m\ = m~ + m 2 q and m 2 _ = m~ — m~ . (3-5) 

3.2 Colour bases in the s-channel 

As discussed in section 2, colour correlations need to be taken into account once NLL 
soft-gluon resummation is performed for processes involving pair-production of coloured 
particles. To this end an appropriate colour basis has to be chosen. We have opted to use 
an s-channel colour basis, which traces the colour flow through the s-channel and has the 
virtue of rendering the anomalous dimension matrices diagonal at threshold [40, 26, 27, 29]. 

Since we are dealing with two coloured particles in both initial and final state, the s- 
channel basis is obtained by performing an s-channel colour decomposition of the reducible 
two-particle product representations into irreducible ones. For squark-squark and squark- 
gluino production this amounts to the following decompositions in SU(3): 

qq^qq: 303 = 306, 

qg^qg: 308 = 3060 15, (3.6) 

where the product representations apply to both the initial and final state. In a general 
SU(iV c )-theory the dimensions of the various representations are of course different, but 
the number of base tensors for these two processes remains the same. 

An economic way to construct the s-channel colour bases for squark-squark and squark- 
gluino production is to start with an arbitrary complete colour basis of the considered 
process in terms of which the s-channel base tensors c/(ai, 02; as, CI4) can be expressed. 
Then the s-channel basis can be obtained by simply requiring that a particular base tensor 
is orthogonal to all other base tensors and projects on itself when contracted in s-channel: 

^2 ci(a 1 ,a 2 ;b, b')cp(b,b'; 0,3,0,4) = ZS n > c/(ai,a 2 ;a 3 ,a 4 ) , (3.7) 

6,6' 

where Z is an arbitrary normalization constant. A similar procedure was found by the 
authors of Ref. [29] on the basis of an analysis in terms of Clebsch-Gordon coefficients. 
This projective construction of the s-channel base tensors constitutes a direct way of ob- 
taining explicit implementations of the irreducible representations on the right-hand side 
of Eq. (3.6). The minimal requirement for the projective method to work is that the par- 
ticles in the initial state must be in the same representations as those in the final state, as 
follows directly from the fact that the labels of the initial state are contracted with those 
of the final state in Eq. (3.7). This is indeed the case for both the squark-squark and 
squark-gluino production processes. An example of the calculation of the s-channel colour 
basis for the qq — > qq process is given in appendix B. 
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In order to present the s-channel base tensors in the subsequent text, we will need 
the following SU(iV c )-objects: the singlet colour structures S a b, where a and b belong 
to particles in either the adjoint or the fundamental representation, the generators of the 
fundamental representation T£ b , the structure constants / a & c and the symmetric forms d a b c - 

3.3 Leading-order partonic cross sections 

Having defined all necessary ingredients, we can now present the results for the colour- 
decomposed qq and qg partonic cross sections at LO. These partonic cross sections are 
averaged over initial-state spin and colour. The colour-decomposed LO cross sections for 
the qq — > qq , gg — ► qq , qq — ► gg and gg — > gg processes, together with their Mellin-moment 
transforms, can be found in Ref. [27]. 

3.3.1 Squark-squark production 

We consider the process 

Qh(ai,Pi)qf 2 (a2,P2) -> q(a3,p 3 ) q{ai,PA) , (3.8) 

where the flavours of the initial-state quarks are indicated by f± , ji and all external particles 
are in the fundamental representation of SU(7V C ). The method described in section 3.2 to 
obtain a suitable s-channel colour basis yields the following two colour tensors: 

C\ = &a\an&a,2az ^0103^0204 and (S^ = 5 aiai 5 a2 a A ~\~ ^0103^0204 ■ (3-9) 

The dimensions of the representations spanned by these two base tensors are given by 
dim(-Rf) = \N C (N C - 1) and dim(i?f ) = \N C (N C + 1). In the SU(3) case this basis 
coincides up to normalization factors with the base tensors given in Ref. [29] for the 3 
and 6 representations. The decomposition of the LO partonic squark-pair cross section in 
terms of the base tensors (3.9) is given by 



(o) 



7raUN 2 -l)(N c + l) 



S 



2m? 2m 2 + s 2mi + sm\ 



2m 2 _ + s /1/2 s " mi + smj 



ircZM - l)(N a - 1) 



W-OT.2 4 N 3 S 



- 2m? 2m 2 + s 2mi + sm? 

2m 2 _+s 1 hh ~ s 1 mt+sml 



with 



' s + 2m^ - sj3 
L | = IO; ' : ' s + 2m 2 _ + s(3 



The quantities (5 and m?_ are defined in Eqs. (2.2) and (3.5), using 1JI3 = 1714 = rriq. The 
occurrence of the Kronecker-delta Sf x f 2 reflects the fact that for equal- flavoured initial-state 
quarks extra diagrams contribute. In appendix A we present results for the Mellin-moment 
transforms of these colour-decomposed LO cross sections. 
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3.3.2 Squark-gluino production 

At the partonic level the qg production process is given by 



Qfi{a>i,Pi)g{a>2,P2) -> 9(03,^3)5(04,^4)- 



(3.10) 



The initial and final state of this process involves both a particle in the fundamental rep- 
resentation (q or q) and a particle in the adjoint representation (g or g). For the s-channel 
colour decomposition the following three base tensors are used: 



J19 



„Q9 



^rpCL4 rpCl2 ~j 



N — 2 

tic ±x x _ or] T c 



N, 



+ 2- c 



c 

c « _ ^C+l 8 X 
c 3 ~ u a 2 a 4 u ai a 3 



4- 2d T c 



C 



N c , 
N c + 2 
N c + 1 



0301 , 



la 3 a\ ■ 



(3-11) 



The dimensions of the representations spanned by these three base tensors are given by 
dim(i?f) = iV c , dim(i?f) = i7V (7V c +l)(iV -2) and dim(i?f ) = ^N C (N C -1)(N C + 2). 
In the SU(3) case this basis coincides up to normalization factors with the base tensors 
given in Ref. [29] for the 3, 6 and 15 representations. The decomposition of the LO 
partonic squark-gluino cross section in terms of the base tensors (3.11) is given by 



a 



(o) 



a 2 ir 



(N 2 - 1)8 



2m?m 2 _ 



2mr + s + 2m? s 9 
— ' — N 2 

2s 2 c 



+ 



ml 



m_ 



sN 2 



+ 



a 



(0) 

Q9^Q9,2 



a (0) 
qg->qg,3 

with 



a 2 MN c - 2) 
(N c - 1)8 

a 2 s Ti{N c + 2) 
(N c + 1)8 



log 



2m\ 
s 

,2 (^2 



7m± +3s Ar2 

A N C 

As c 



L 2 

3m 2 , + s 



2s 



+ 



7m 2 



4N 2 s 



k(5 



2m_(m 2 \_ — s) — s 2 
Is~ 2 

2m 2 _{m 2 + - s) - s 2 
I? 

s + m 2 . — ks(3 



m (mi — s) m 
L2+ 2s* L *~—^ 



L 2 + 



2 / 2 

m_(m. 



s + m_ + Ksj3 



and 



La 



2s 2 



log 



s m_ 

— L 3 k/3 

s 



s — m_ 



Ksfi 



s — m_ + ks(5 



The quantities (5 , k and m± are defined in Eqs. (2.2), (3.4) and (3.5), using 7713 
and 711,4 = m-g. In appendix A we present results for the Mellin- moment transforms of 
these colour-decomposed LO cross sections. 



3.4 The soft anomalous-dimension matrices 

As we reviewed in section 2 below Eq. (2.11), resummation to NLL accuracy requires the 
anomalous dimensions Tjj of the products of Wilson-line operators connected by a base 
tensor c/. To this end one must compute the UV divergences from their loop corrections, 
and from these the renormalization constants Zjj for these operators. Here we only need 
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the one-loop corrections. The anomalous dimensions can be computed from the residues 
of the UV poles in the renormalization constants Zjj as 

d 

Tu = -a s — Res e ^o^/j(a s ,e). (3.12) 

The relevant UV divergences occur in loop corrections to the base tensors cj [40, 37] due 
to the Wilson lines. The complete first order correction to cj can be written as 

Y,u ij C?jCj, (3.13) 

ij 

where i and j denote the eikonal lines between which the gluon is spanned, uj l i is the 
corresponding kinematic part of the one-loop correction, and C l j j denotes how the base 
tensors get mixed due to the corrections. At one-loop we can calculate the anomalous 
dimensions directly from Eq. (3.13) 

The precise form of this function depends on the colour basis chosen. The eikonal integrals 
that constitute the uj ij can be found in Ref. [40], except for the unequal-mass case that 
we need for squark-gluino production. The corresponding integral uj 34 is discussed in 
appendix D, using the Feynman rules in the eikonal approximation presented in appendix C. 

In order to present the results for the soft anomalous dimensions in a compact way, 
we introduce the following t- and ^-channel quantities 

A = i [T(m 3 ) + T(m 4 ) + U(m 3 ) + U(m 4 )] , 

n = l[T(m 3 ) + T(m 4 )-U(m 3 )-U(m 4 )], (3.15) 

in terms of the t- and u-channel logarithms 3 

m . . , / m 2 — 1\ 1 — m , . . , ( m 2 — u\ 1 — in 

T{m) = log -== — and U(m) = log — . (3.16) 

\ vsm z J z \ vsm z J £ 

The one-loop soft anomalous-dimension matrices for the qq — > qq , gg — > qq , qq — > gg 
and gg — ► gg processes have been calculated in Ref. [27], where the corresponding values 
of the Dij^ki j coefficients can be found as well. 4 

3.4.1 Soft anomalous dimensions for squark-pair production at one- loop 

In the basis (3.9) the one-loop soft anomalous-dimension matrix is given by 



f 2l 



(c 2 {Rf)K- ^£±1(L P + 1) -(N c + l)n \ 



^ -(N c -l)n c 2 {Rf)A + \-±{L p + l) j 



, (3.17) 



3 Note that in the case of equal masses m-j, = 7714 the quantities A,Q,,T(m) and U(m) reduce to the 
corresponding quantities A, Cl,T and U defined in Ref. [27]. 

4 Note that Ref. [27] uses a subtraction term different from Eq. (2.11). 
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with 



2(3 



log([ T |)+,vr 



The coefficients C 2 {R qq ) for I = 1,2 are the quadratic Casimir invariants belonging to 
the representations spanned by the base tensors c qq : 

C 2{ Rf) = (^ + l)(^-2) and C2{Rf) = (N c -l)(N c+ 2)_ (318) 



3.4.2 Soft anomalous dimensions for squark-gluino production at one-loop 

In the basis (3.11) the one- loop soft anomalous-dimension matrix is given by 

( ^ 



qg^qg 



2vr 



3« 



4iV2(iV c -2) 



ft 



4iV2(iV c + 2) 



(7V2 _ l)(iV - 1) (AT2_i )(iVc + l) 



N c (N c -2) 
iV c -l 







N c + 1 

f 33, q 9 



ft 



(3.19) 



with 
Til, <? 9 

^22,99 
r 33, 99 

where 



1 iV 2 + 1 

C 2 (i?f ) A + [C F + — ] n - [ T ( m .) _ r(m -)] _ at c (L, 3i „ 4 + 1) , 



2N, 



C 2 {Rf)A + [C F 
C 2 (Rf)A + [C F 



N, 



— ]n 



1 



JV„ + 1- 



ft 



c 

iVg + 1 
2iV c 

^c + ! 

2JV„ 



[T(m 9 -) - T(m 5 )] - (L V3>V4 + 1) , 



[T(mg) - T{m~g)] + (^ 3jl , 4 + 1) , 

(3.20) 



J V 3 ,VA 



2k(3 



. k — (3 



(3.21) 



The explicit derivation of Eq. (3.21) is presented in appendix D. The coefficients C 2 (R q9 ) 
for / = 1,2,3 are the quadratic Casimir invariants belonging to the representations 
spanned by the base tensors c q9 : 



C 2 (R q 



qg\ 



2JV„ 



= C F , C 2 (R q9 ) 



(N c -l)(3N c + l) 
2N^ 



and C 2 (Rf) 



(JV + 1)(3JV C -1) 
2AL 



(3.22) 
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3.4.3 The threshold limit 

At the production threshold, where (3 — > 0, the soft anomalous-dimension matrices become 
diagonal by virtue of using an s-channel basis. In addition, the diagonal components 
become proportional to the total colour charge of the heavy-particle pair produced at 
threshold: 

Dij^kij = -C 2 {Rf) , (3.23) 

with C2(R l j) as given in equation (3.18) for squark-pair production and in equation (3.22) 
for squark-gluino production. In the SU(3) case the Aj-+fc«,/ coefficients for squark-pair 
production are given by 

{Dqg^qgj} = {-4/3,-10/3}, 

while for the squark-gluino production process they are 

{D qg ^q § ,i} = {-4/3,-10/3,-16/3}. 

4. Numerical results 

In this section we present numerical results for the NLL-resummed cross sections matched 
with the complete NLO results for squark and gluino pair-production at both the Tevatron 
(y/S = 1.96 TeV) and the LHC (y/S = 14 TeV). The matching is performed according to 
Eq. (2.13). From now on we refer to the matched cross sections as NLL+NLO cross sections. 
We also compare the NLL+NLO predictions with the corresponding NLO results. The NLO 
cross sections are calculated using the publicly available PR0SPIN0 code [44], based on the 
calculations presented in Refs. [14, 15, 16]. As described in detail in Ref. [16], the QCD 
coupling a s and the parton distribution functions at NLO are defined in the MS scheme 
with five active flavours. The masses of squarks and gluinos are renormalized in the on- 
shell scheme, and the SUSY particles are decoupled from the running of a s and the parton 
distribution functions. As already discussed in previous sections, no top-squark final states 
are considered. We sum over squarks with both chiralities (qi and q~ii), which are taken as 
mass degenerate, and include the charge-conjugated processes in the numerical predictions. 
For convenience we define the average mass of the sparticle pair m = (m^ + m4)/2, which 
reduces to the squark and gluino mass for qq, qq and gg final states, respectively. The 
renormalization and factorization scales fx are taken to be equal. In order to evaluate 
hadronic cross sections we use the 2008 NLO MSTW parton distribution functions [43] 
with the corresponding a s (M§) = 0.120. The numerical results have been obtained with 
two independent computer codes. 

We first discuss the scale dependence of the NLL+NLO matched cross section for the 
separate processes pp — ► qq , qq , qg , gg + X at the Tevatron. Figure 1 shows the NLO and 
NLL+NLO cross sections for niq = m g = m = 500 GeV as a function of the renormalization 
and factorization scale fi. The value of (i is varied around the central scale no = m from 
[i = no/ 10 up to [i = 5 fiQ. As anticipated, we observe a reduction of the scale dependence 
when going from NLO to NLL+NLO, in particular for gg and qg production (Figs, lb 
and Id, respectively). In the case of squark pair-production, on the other hand, the scale 
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reduction due to soft-gluon resummation is moderate (see Figs, la and lc). We note that 
the gluino-pair production cross section (Fig. lb) is rather small for this particular choice 
of masses because of a suppression of the LO qq — ► gg amplitude proportional to mj — m| 
near threshold (cf. Eq. (55) of Ref. [16]). 

At the central scale fi = fiQ = m the cross-section predictions are in general enhanced 
by soft-gluon resummation. The relative if-factor Knll — 1 = 0"nll+nlo/Vnlo — 1 at 
the Tevatron is displayed in Fig. 2 for squark and gluino masses in the range between 
200 GeV and 600 GeV. We show results for various mass ratios r = nig/nig. The soft-gluon 
corrections are moderate for qq production (Fig. 2a), but reach values up to 27%, 29% and 
60% for gg , qq and qg final states, respectively, in the range of r we consider. Because of the 
increasing importance of the threshold region, the corrections in general become larger for 
increasing sparticle masses. The strong r-dependence of -K'nll for gluino-pair production 
in Fig. 2b is driven by the r-dependence of the NLO cross sections for qq — ► gg. The 
large effect of soft-gluon resummation for qg and gg production can be mostly attributed 
to the importance of gluon initial states for these processes. Furthermore, the presence 
of gluinos in the final state results in enhancement of the NLL contributions [27], since in 
this case the Casimir invariants that enter Eq. (2.8) reach higher values than for processes 
involving only squarks. The substantial value of i^NLL fo r QQ production at the Tevatron 
is a consequence of the behaviour of the corresponding NLO corrections, which strongly 
decrease with increasing squark mass [16]. 

We now turn to the discussion of pair production of squarks and gluinos at the LHC, 
i.e. pp — ► qq , qq , qg , gg + X. The results for the processes pp — > qq and pp — > gg agree with 
those presented in Refs. [26, 27], while the predictions for pp — > qq and pp — ► qg are new. In 
Fig. 3 the cross sections are shown for squark and gluino masses nig = m g = m = 1 TeV as 
a function of the common renormalization and factorization scale [i. The scale uncertainty 
of the theoretical prediction is reduced at NLL+NLO. Similarly to the Tevatron case, 
soft-gluon resummation is most significant for gluino-pair production and squark-gluino 
production. For those processes, the relative if -factor Knll — 1 reaches 35% for gluino- 
pair production and 18% for squark-gluino production at the highest accessible sparticle 
masses around 3 TeV (see Figs. 4b and 4d). The r-dependence of i^NLL f° r gluino-pair 
production is again driven by the r-dependence of the NLO cross section, discussed in 
Ref. [16]. 

Representative values for the NLO and NLL+NLO cross sections at the Tevatron and 
the LHC are collected in Tables 1 and 2 for equal squark and gluino masses. 

The impact of the NLL resummation on the cross section for inclusive squark and 
gluino production, i.e. pp/pp — > qq + qq + qg + gg + X , can be inferred from the inclusive 
-factor displayed in Fig. 5. The pattern exhibited in Fig. 5 can be understood from 
the relative importance of the qq ,qq , qg and gg final states and from their individual 
if-factors as shown in Figs. 2 and 4. At rriq = rrig « 400 GeV, for example, the inclusive 
cross section at the Tevatron (Fig. 5a) is built up from the individual final states in the 
ratio qq : gg : qg : qq ~ 1 : 3.6 : 14 : 32, as can be read off from Table 1. Owing to 
the large NLL corrections for the qg final state, the resulting inclusive if -factor JTnll is 
approximately 1.1. At rriq = m g = 600 GeV the correction to the inclusive cross section at 
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the Tevatron due to NLL resummation can be as high as 18%. The inclusive corrections are 
smaller at the LHC for sparticle masses below 3 TeV (see Fig. 5b). Given the sparticle mass 
ranges that we consider, this is consistent with the fact that the distance from threshold, 
i.e. the value of the variable 1-/9 = 1 — 4m 2 /S, is on average larger at the LHC than at 
the Tevatron. 

In Figs. 6a and 6b we show for the Tevatron and LHC, respectively, the resummed 
NLL+NLO total cross section for inclusive squark and gluino production as a function of 
the average sparticle mass m. For illustration we show these results for the choice rriq = m g . 
The error bands indicate the theoretical uncertainty of the NLL+NLO total cross section 
due to the scale variation in the range m/2 < fi < 2m. The results presented in Fig. 6 
are the most accurate theoretical predictions currently available for the above processes. 
The reduction of the theoretical error due to variation of the common factorization and 
renormalization scale /x between /i = m/2 and fi = 2m is illustrated in Fig. 7a for the 
Tevatron and in Fig. 7b for the LHC. Both at the Tevatron and at the LHC, soft-gluon 
resummation leads to a significant reduction in this part of the theoretical uncertainty. 



5. Conclusions 

We have performed the NLL resummation of soft gluon emission for squark and gluino 
hadroproduction. Explicit analytical results are presented for the anomalous dimension 
matrices and the colour-decomposed LO cross sections in x and iV-space for the qq and qg 
final states. We provide NLO+NLL matched numerical predictions for all pair-production 
processes of coloured sparticles at the Tevatron and the LHC. The NLL corrections lead 
to a significant reduction of the scale dependence and, in general, increase the NLO cross 
sections. The effect of soft-gluon resummation is most pronounced for processes with 
initial-state gluons and final-state gluinos, which involve a large colour charge. Specifi- 
cally, at the Tevatron we find an increase of the cross-section prediction of up to 40% at 
sparticle masses around 500 GeV when going from NLO to NLL+NLO, depending in de- 
tail on the final state and the ratio of squark to gluino masses. For the inclusive sparticle 
cross section at the Tevatron, summed over all pair-production processes for squarks and 
gluinos, the enhancement can be as large as approximately 15% in the mass range up to 
500 GeV, probed by current experimental searches. At the LHC, the NLL corrections are 
particularly significant for squark-gluino production and gluino-pair production, reaching 
approximately 20% and 30%, respectively, for sparticle masses around 3 TeV. Both at the 
Tevatron and at the LHC, the inclusion of NLL corrections leads to a reduction of the scale 
dependence over the full mass range that will be probed by experiments. In addition, the 
NLL corrections lead to a significant enhancement of the NLO cross-section predictions for 
heavy sparticles. The NLL+NLO matched predictions presented in this paper should thus 
be used to interpret current and future searches for supersymmetry at the Tevatron and 
the LHC. 
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Figure 1: The scale dependence of the NLL+NLO and the NLO total cross sections for squark 
and gluino pair-production processes at the Tevatron. The squark and gluino masses have been set 
to rriq = m~g = m = 500 GeV. 
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Figure 2: The relative NLL if-factor Knll — 1 = onll+nlo /^nlo — 1 f° r squark and gluino 
pair-production processes at the Tevatron as a function of the average sparticle mass m. Shown are 
results for various mass ratios r = mg/rriq. 
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Figure 3: The scale dependence of the NLL+NLO and the NLO total cross sections for squark 
and gluino pair-production processes at the LHC. The squark and gluino masses have been set to 
rriq — nig — m — 1 TeV. 
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Figure 4: The relative NLL if-factor -STnll — 1 = ctnll+nlo /cnlo — 1 f° r squark and gluino 
pair-production processes at the LHC as a function of the average sparticle mass m. Shown are 
results for various mass ratios r = mg/rriq. 
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Figure 5: The relative NLL AT- factor ATnll — 1 = ctnll+nlo/cnlo — 1 f° r the inclusive squark 
and gluino pair-production cross section, pp/pp — > qq + qq + qg + gg + X, at the Tevatron (a) and 
the LHC (b) as a function of the average sparticle mass m. Shown are results for various mass 
ratios r — mg/rriq. 
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Figure 6: The NLL+NLO cross section for inclusive squark and gluino pair-production, pp/pp — > 
qq + qq + qg + gg + X , at the Tevatron (a) and the LHC (b) as a function of the average sparticle 
mass to. Shown are results for the mass ratio r — mg/mq = 1. The error band corresponds to a 
variation of the common renormalization and factorization scale in the range to/2 < /x < 2m. 
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Figure 7: Scale dependence of the NLL+NLO and NLO cross sections for inclusive squark and 
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Table 1: The NLL+NLO and NLO cross sections for the squark and gluino pair-production pro- 
cesses at the Tevatron. Shown are results for the mass ratio r = mg/m q = 1. The common 
rcnormalization and factorization scale has been set to to. 
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1.14 x 10" 


-3 


7.59 x 10" 


-6 


Knll - 1 


0.010 


0.012 


0.017 


0.034 




0.064 






pp -> gg at 


^S= 14 TeV (r=1.0) 










mg [GeV] 


200 


500 


1000 


2000 


3000 


cnlo [pb] 


3.74 x 10 3 


2.85 x 10 1 


2.92 x lO" 1 


5.82 x 10" 


-4 


2.68 x 10" 


-6 


cnll+nlo [pb] 


3.86 x 10 3 


3.00 x 10 1 


3.18 x 10" 1 


6.91 x 10" 


-4 


3.62 x 10" 


-6 


Knll — 1 


0.033 


0.054 


0.089 


0.19 




0.35 






pp -> qq at y/S = 14 TeV (r=1.0) 










rriq [GeV] 


200 


500 


1000 


2000 


3000 


ctnlo [pb] 


5.45 x 10 2 


1.34 x 10 1 


5.28 x 10" 1 


6.48 x 10" 


-3 


1.18 x 10" 


-4 


cnll+nlo [pb] 


5.46 x 10 2 


1.34 x 10 1 


5.32 x 10" 1 


6.64 x 10" 


-3 


1.25 x 10" 


-4 


Knll — 1 


0.003 


0.004 


0.008 


0.024 




0.056 






pp -> ^ at 


s/S = 14 TeV (r=1.0) 










to [GeV] 


200 


500 


1000 


2000 


3000 


TNLO [pb] 


4.86 x 10 3 


6.55 x 10 1 


1.22 


5.49 x 10" 


-3 


4.96 x 10" 


-5 


cnll+nlo [pb] 


4.92 x 10 3 


6.69 x 10 1 


1.26 


5.96 x 10" 


-3 


5.80 x 10" 


-5 


Knll — 1 


0.013 


0.021 


0.037 


0.085 


0.17 



Table 2: The NLL+NLO and NLO cross sections for the squark and gluino pair-production pro- 
cesses at the LHC. Shown are results for the mass ratio r = TOg/m, = 1. The common renormal- 
ization and factorization scale has been set to to. 
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A. Leading-order AT-space cross sections for qq and qg production 



In this appendix we present the analytical results for the Mellin transforms of the LO 
cross sections for qq and qg production. The cross sections are colour-decomposed in 
SU(3) according to the procedure described in section 3. The Mellin-transformed LO cross 
sections for the qq and gg final states can be found in [27]. 

The expressions for the colour-decomposed LO A^-space cross sections for the process 
9/i 9/ 2 99 are given by 



*(°) 



qq-^qq 



27m? 



- 5 flf2 H N - 



4B N G N 
2N + 3 



f 2r 2 1 \ 

V + r 2 + 1 TV + 2 ) 



+ 2B N 



N 2 + 2N + 2 
{N + l)(N + 2) 



27m? 



2B n Gn 



[N + -= 

V r 2 + 1 N + 2 J 



(A.l) 



+ B 



N 



N 2 + 2N + 2 
{N + l)(N + 2) 



(A.2) 



whereas for the process qg — > qg they read 



"W^WAVV 8m 2 



9B N+1 P N+l (l-r) _ 9B N P^ B N+2 P N+2 (7r 2 - 9)(1 
.•T- - - , ,, + (r + 1)5 



-r) 



(r + 1)3 



2(r + l) 2 



+ 



+ 



Bn+iP£ +1 (1 ~ r) _ B N+2 P+ +2 (r 2 + 17)(l- 
9fr + l) 3 9(r + l) 5 



130B N +1 K N+1 (l-r) _ 56B N K N 
1 " :! 9(r + l) 2 



2 

8m] 



9(r + 1) 

2B N+1 P N+1 (l-r) B N P N 2B N+2 P N+2 (r 2 + l)(l 
(r + l)3 (r + 1) 2 (r + 1) 5 



(A.3) 
-r) 



+ 



+ 



2£? JV+1 P+ +1 (l-r) _ 2B N+2 P+ +2 (r 2 + l)(l-r) 
(r + 1) 3 (r + 1) 5 

4B N+1 K N+l (1 - r) 
(r + 1) 3 



(A.4) 



(A.5) 



- 25 - 



We have used the following abbreviations: 
B N = p(N + 1,1/2), 



G N = 2 Fi 1,1/2,^ + 5/2, 



2 1 \ 2N 

r z — 1 



r 2 + 1 



P± = _1- 2 F X ( l/2,iV + l,iV + 3/2, ' "' 



x ~ N + 1 
1 -r 



r + lj N + 3/2 
1 



1 + r 



' 2^(1/2,^ + 2,^ + 5/2/^ 



2iV + 3 



2 Fx -l/2,JV + l,JV + 5/2, 



1 — r 



if 



A' 



dz- 



1 / l-r 



2r- 



log 



1 + r 

2(1 + y 7 !"^) + (r 2 - l)z 
2(1 - v 7 !^) + (r 2 - l)z 



(A.6) 



with 2 -Pi (A, fj,, u, £) the hyper geometric function, f3(fi, u) the beta function and r = rrig/mq. 
For the numerical evaluation of Hn we use the expansion 

(-l) k 



9t ,2 / r 2 1 \ '« 1 

^ = -^W^4) — y 



m=0 

13 (k + N + 2, 1/2) 
/c + iV + 2 

2 iq [ 1,1/2,* + AT + 7/2, 



r 2 - 1 
1 + r 2 



l-r' 
1 + r s 



/?(& + l,m- k + 1) 



/3(£; + iV + 2,3/2) 



2\ 2^ 



(A.7) 



B. Construction of the s-channel colour basis: an example 




Figure 8: The LO diagrams that contribute to squark-pair production. 



For the process qq — * qq we explicitly show how to derive the s-channel colour basis 
given in Eq. (3.9). The same steps can be used to obtain the basis given in Eq. (3.11) for 
the qg — ► qg process, although the calculations are more tedious in that case. 
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As a starting point we take the colour structures that occur in the LO qq — > qq process 
displayed in Fig. 8. Using the conventions introduced in section 3.1 these are: 



rpc rpc 
£13(11 (24(22 



rpc rpc 
(23(22 (24(21 



2 ^0302^04(21 ^0301^0402 ] ' 

2 ( ^0301^0402 ~~ ^0302^0402 1 ■ 



where c is a summation index in the adjoint representation. For convenience these colour 
structures have been rewritten in terms of the t and it-channel singlet structures S a3ai 5 a4a2 
and 5 a3a2 5 aiai . It is clear from this expression that two independent singlet structures 
occur. Since the two-particle reducible product representation 3®3 contains two irreducible 
representations, cf. Eq. (3.6), this basis must be complete. That means that the s-channel 
base tensors are linear combinations of these singlet structures. The projective prescription 
(3.7) leads to the following set of equations: 



^03 b' "04 b + Bi'5a 3 bS ai b') = Z5ir(AiS 

0302^0401 Bj5 a3ai 5 aia2 




where 1,1' £ {1,2} and Z is an arbitrary normalization constant. Working out the equa- 
tions shows that up to interchanging the base tensors the unique solution is given by 
A\ = —A2 = B% = £>2 = Z/2, which is exactly the basis given in Eq. (3.9). 

One can check explicitly that this basis is complete 
for gluon resummation: representing the combined colour 
structure of the external particles by one of the base ten- 
sors cr and connecting any two external particles by an 
additional gluon yields no additional colour structures. In 
Fig. 9 an example of such a gluon insertion is shown. For 
processes for which the LO colour basis is not complete, 
this procedure can also be used to identify additional base 
tensors. 

If a particle is exchanged in the s-channel, the cor- 
responding base tensor has a direct physical interpre- 
tation. An example is the Feynman diagram for the 
qg — * qg process shown in Fig. 10. Since the quark 
exchanged in the s-channel is in the fundamental repre- 
sentation, the corresponding N c -dimensional base ten- 
sor (cf in Eq. (3.11)) can be read off immediately from Figure 10: Example of a diagram 
the colour structure of this diagram. corresponding to a base tensor. 



Figure 9: 

insertion. 



An example of gluon 




C. Eikonal Feynman rules 

In this appendix the eikonal Feynman rules will be given for a soft gluon with momentum 
k attached to an eikonal line with momentum p. In the eikonal approximation we have 
k <C p, which leads to simple Feynman rules since the propagator that connects the matrix 
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element to the radiated gluon becomes effectively on-shell. The generic diagrams and their 
corresponding Feynman rules are given by (cf. [40]) 



p p 



° = 9s{TCRU p-irr- -^O (ci) 



[1,C 

for an incoming eikonal line and 



J) 



p p 



O = 9s(T c R ) ab -i^— Ql (C.2) 




p ■ k + ie 



for an outgoing eikonal line. Here g s is the strong coupling constant, \x is the Lorentz 
index of the gluon and ie represents the infinitesimal imaginary part of the propagator 
that connects the matrix element to the radiated gluon. The colour labels of the different 
particles are denoted by a, b and c. The representation of the eikonal line is denoted by 
R. We have R = F for the fundamental representation, R = F for the charge conjugate 
of the fundamental representation, and R = A for the adjoint representation. The colour 
operators occurring in Eqs. (C.l) and (C.2) are given in Table 3. Note that the order of 
the colour indices a, b, c in f a b c is kept fixed irrespective of whether the gluon is emitted 
above or below the eikonal line. 

Outgoing (s)quark / incoming anti-(s)quark: (Tp) a b = T^ b 
Outgoing anti-(s)quark / incoming (s)quark: (Tp) a {, = —T^ a = —(T^)* 
Gluons / gluinos: {T%)ab = Kb = -ifaJbc 

Table 3: Colour operators used in the eikonal Feynman rules. 



D. One-loop eikonal integral for qg production 

We briefly present here the calculation of the kinematic part w 34 of the one-loop correction 
to the process qg — > qg in the eikonal approximation. The equal-mass case of lo is well 
known [40], but for qg final states we also need the unequal-mass version. 

The kinematic part of the one-loop correction generated by the exchange of a virtual 
gluon between the two final-state eikonal lines is according to Eq. (C.2) given by 

- 34 = 9lfP^( ) ■ ( ~ h ) ttP-, N^(k) . (D.l) 

Hs J (2ir) d \v 3 -k + ieJ \-v 4 -k + ieJ (k 2 + ie) y 1 y ' 

We use dimensionless vectors v 1 ^ = W2/s with pi denoting the momentum of the massive 
external particle i. We calculate the gluon propagator in a general axial gauge with 

n ■ k (n ■ k) z 
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where is a general gauge vector with n 2 < 0. In the case that v 2 4 > and u| / v\ the 
solution of the integral w 34 reads 



LO 



34 



Q: 



7re 



V4 



1] , 



(D.3) 



with e = 4 — d. The gauge-independent term L VZ)VA is given by 



L 



V 3 ■ V 4 



V3,V4 



2 (v 3 • v A ) 2 - v\v\ 



2m + log 



2„,2 



3^4 



+ V3 ■ f 4 + V 7 ( w 3 • V4) 2 - V3V4 



+ log 



•^3 + V3 ■ VA ~ \J (V3 ■ V4) 2 ~ vjvj 

vl + v 3 ■ v 4 + v 7 i v 3 • V4) 2 ~ v 2 vj 



(D.4) 



The gauge-dependent terms L V3 and L^ 4 can be found in Ref. [40] and cancel against 
contributions from the self-energy diagrams when calculating the anomalous dimensions. 
The gauge- independent term L V3>m can be rewritten in a compact form using and k as 
defined in Eqs. (2.2) and (3.4): 



L 



V3,V4 



K 2 +f3 2 

2k(3 



+ in 



(D.5) 



For equal- mass final-state particles this quantity reduces to the well-known form (cf. [27]) 



La = 



l + (3 2 
213 



log 



+ in 
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